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Abstract 

The problem of estimating the parameters of a linear regression model Z(s, t) = 

Wi5i(s, t) H \-Tn p g p (s, t) + U(s,t) based on observations of Z on a spatial domain 

G of special shape is considered, where the driving process U is a Gaussian random 
field and gi, . . . ,g p are known functions. Explicit forms of the maximum-likelihood 
estimators of the parameters are derived in the cases when U is either a Wiener 
or a stationary or nonstationary Ornstein-Uhlenbeck sheet. Simulation results are 
also presented, where the driving random sheets are simulated with the help of their 
Karhunen-Loeve expansions. 
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1 Introduction 



The Wiener sheet is one of the most important examples of Gaussian random fields. It has 
various applications in statistical modeling. Wiener sheet a ppears as limiting; process of some 
random fields defined on the int erface of the Is ing model (IKuroda and Manakal . 119871 ). it is 
used to model random polymers (iDouglasl . I1996T). to describe the dynamics of Heath- Jarrow- 



Morton t ype forward interest ra t e mod els (IGoldsteinl. 12000) or t o model random mortality 



surfaces (jBiffis and Millossovichl . 120061 ). Further, ICarterl (120061 ) considers the problem of 
estimation of the mean in a nonparametric regression on a two-dimensional regular grid of 
design points and constructs a Wiener sheet process on the unit square with a drift that is 
almost the mean function in the nonparametric regression. 
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The stationary Ornstein-Uhlenbeck sheet {X(s,t) : s,t 6 R} is a zero mean Gaussian 
process with covariance structure 

J2 



EX(s l ,t 1 )X(s 2 ,t 2 ) 



a 



4a(3 



-ct|s 2 — si|— /3|t 2 — *l| 



'1.1) 



where a > 0, /3 > 0, a > 0, while the random field 

s t 





;i. 2 ) 



where a6l, G R, cr>0 and {W(s, £) : s,t > 0} is a standard Wiener sheet, can be 
considered as the Ornstein-Uhlenbeck sheet with zero initial condition on the axes. 



Ornstein-Uhlenbeck sheets play role e.g. in potential theory (IFeyel and de La Pradelle 



19951 ) and, si milarly to the Wiener sheet, they also appea r as dr iving fields in forward interest 
rate models (jGoldsteinl . 12000c ISanta-Clara and Sornettd . 1200 ll ). 

In this paper we consider a linear regression model driven by a Gaussian sheet, that is a 
random field 

Z(s, t) := migi(s, t) H h m p g p (s, t) + U (s, t) (1.3) 

observed on a domain G, where gi,...,g p are known functions and U is either a Wiener or 
an Ornstein-Uhlenbeck sheet, and we determine the maximum likelihood estimator (MLE) 
of the unknown parameters mi, . . . , m p . 

In principle, the Radon-Nikody m derivative of Gaussian measures might be derived from 
the general Feldman-Hajek theorem (lKuol . ll975l ). but in most of the cases explicit calculations 
can not be carried out. In the case when U is a standard Wiener sheet, p — 1 an d qi = 1 



Rozanov 



(shifte d Wiener sheet), the MLE of the unknown parameter is given e.g. in 
(119681 ). where the estimator is expressed as a function of a usually unknown random variable 
satisfying some characterizing equation. In several cases the exact form of this random 



variable can be derived by a method proposed bylRozanovi (1l990l ). based on linear stochastic 
partial differential equations. lArato. N. M.I (Il997al ) used Rozanov's method to find the MLE 



of the shift parameter of a shift ed Wiener she e t obser ved on a special domain. iBaran et al. 



(120041 ) considered the model of lArato. N. M.I (Il997al ). and applying an essentially simpler 



direct discrete approximation approach the authors found the MLE of the shift parameter 
under much weaker conditions. Later this discrete approximation was used to determine the 
MLE of the unknown parameter for the model (11.31) with p = 1 and a more complicated 
domain of observations, when U is a standard Wiener sheet. 



Arato. N. M.I (jl997bl ) also studied the case when U is an Ornstein-Uhlenbeck sheet, 
p = 1 and g\ = 1, and using partial stochastic differential equations found the MLE of 
the unknown parameter based on the observation o f the r andom field Z on a rectangular 
domain. This result was generalized by lBaran et al\ (120031 ) for the case p—1 and g± with 
slight analytic restrictions. 
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In the present paper we consider the same type of domain G as in iBaran et all (120111 ) 
and extend their result for the general model (11. 3p . We also consider the cases when the 
driving proces s U is a stati o nary a nd a zero start Ornstein- Uhlenbeck sheet and generalize 
the results of lArato. N. M.I ( 1l997bl ) and IBaran et al\ (120111 ). Moreover, we present some 
simulation results to illustrate the theoretical ones where the driving Gaussian random sheets 
are simulated with the help of their Karhunen-Loeve expansions. The proofs of the theorems 
are given in the Appendix. 



2 Models and estimators 

Consider the model (II. 3ft with some given functions g±, . . . , g p : R+ — > R and with unknown 
regression parameters mi, . . . ,m p G R. Let [a, c] C (0, oo) and bi,b 2 G (a, c), let 
7i,2 : K &i] — > R and 7 : [62, e] — > R be continuous, strictly decreasing functions and 
let 71 : [b\,c] —> R and 72 : [a, 6 2 ] — R be continuous, strictly increasing functions with 
71,2(61) = 11(h) > 0, 72(62) = 70(62), 71,2(0) = 72(a) and 71(c) = 70(c). Consider the 
curve T := i\ 2 U Ti U T 2 U r , where 

Ti.2 := |(s,7i,2(s)) : s G [a, ft]]}, := {(s,7i(s)) : s G [6i,c]|, 

T 2 := |(s,7 2 (s)) : s G [a,b 2 ]\, T := |(s,7 (s)) : s G [6 2 ,c]|, 

and for a given e > let T\ 2 , T\, Y £ 2 and Tq denote the inner e-strip of r 12 , r 1; T 2 
and r , respectively, that is e.g. 

r i,2 : = {(*,t) G R 2 :s G [a,a + e], t G [71,2(5), 71.2(a)] or 

s G [a + e,6i], t G [71,2(5), 71,2(5) +e]}. 

Suppose that there exists an e > such that 

qnr| = and rf 2 nr^ = 0, (2.1) 

and consider the set G := G\ U G 2 U G3, where 

Gx := {(s,t) G R 2 : s G [a, 6 X A 63], t G [71,2(5), 72(5)]}, 
G 2 :-- 



{{s,t) G R 2 : s G [61,62], t G [7l(5), 72(s)]}, if 61 < 62 
{(s,t) e R 2 : s G [62,^], t G [71,2(5), 7o(5)]}, if 61 > 6 2 
G 3 := {(5,t) GR 2 : s G [61 V6 2 ,c], t € [71(5), 70(5)]}. 



An example of such a set of observations can be seen of Figure [TJ 

First we study the case when U a standard Wiener sheet and consider the random 
field Z(s, t) := m 1 (? 1 (s, t) H — ■ + m p g p (s, t) + W(s, t). The following theorem is an obvious 
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a bj b 2 c s 

Figure 1: An example of a set of observations G. 



extension of Theorem 2.1 of iBaran et all ( 120111 ) and can be proved exact ly in the same way . 



The proof is based on the discrete appro ximation met hod described in IBaran et al\ ( 12003 



20041 l201lh . which relies on the results of lAratol (I1982L Section 2.3.2). 



Theorem 2.1 If gx, . . . ,g p are twice continuously differentiable inside G and the partial 
derivatives d±gi, dig% and did 2 gi, i = l,...,p, can be continuously extended to G then 
the probability measures Pz and Pw, generated on C(G) by the sheets Z and W , 
respectively, are equivalent and the Radon- Nikodym derivative of Pz with respect to Pw 
equals 

dP z 



dP 



exp I — — (m T Am — 2( T m) j- 



where A := (A k/ ) P k e=v m := (mi, . . . , m p ) T and ( := (d, . . . , ( P ) T with 
^(61,71,2(61)) ^(61,71,2(61)) 



A 



kl 



In 



+ 



6171,2(61) 

[#fc(s, 7i,2(s)) - sdig k (s,-/ li2 (s))] [^(s,7i )2 (s)) - s<9i<^( s , 7i,2(«))] 



(2.2) 



s 2 7i,2(s) 



ds 



<9i#fc(s,7i(s)) dig e (s,jx(s)) 



72(62) 



ds 



In 



7H«J 



72(a) 



d 2 g k (%\t),t)d 2 g e (^ 1 (t),t) 



dt 



7l,2(a) 



71,2(01) 



digk (71,2 (t) ; *) d 29£ (71,2 (*) ; *) 



dt+ d 1 d 2 gk(s,t)did 2 ge(s,t)dsdt, 



G 
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and 



0fc(&i,7i,2(&i))Z(&i, 71,2(61)) / <9i# fc (s,7i(s)) , , 



6171,2(61) 



In 



bi 



+ 



a 

72(62) 



[gkjs, 71.2(g)) - s^ifffc (a, 71.2(g))] 

g 2 7i,2(s) 
«9 2 ^(7 2 _1 W,0 



7i(g) 

^(s, 71.2(g)) ds - sZ(ds,7i j2 (s)) 

71.2(a) 



(2-3) 



7 2 - 1 W 



z( 72 - 1 (t),dt) 



71,2(01) 



7r 2 i (t) 



z( 7 rj(*),d*) 



72(a) 

+ ^^(S, t)Z(ds,dt). 

G 

T/ie maximum likelihood estimator of the parameter vector m based on the observations 
{Z(s,t) : (s,t) G G} has the form rh = A~ 1 ( and has a p- dimensional normal distribution 
with mean m and covariance matrix A^ 1 . 

Remark 2.2 We remark that the weighted L 2 -Riemann integrals of partial derivatives 
of the Wien e r shee t (and of other L 2 -processes) along a curve are defined in the sense of 
Baran et al\ ( 1201 ll . Definition 4.1). This means that if Z is an L 2 -process given along an 
^-neighborhood of a curve T := {(5,7(5)) : s G [a, 6]}, where 7 : [a, b] — > R is strictly 
monotone and y : [a, b] — y M. is a function, then 

b b 
J y(s) Z{ds, 7(g)) := LLm. \ j y(s) [Z(s + h, j(s)) - Z(s, 7(g))] ds, 



7 (b) 



7 (6) 



y( 7 - 1 (0)Z( 7 - 1 (t),dO :=l.^m.^ / v(7 _1 (*)) P(7 _1 (*),* + *) ~ ^(7 _1 (*).*)] d*. 



7(a) 



7(a) 



if the right hand sides exist. 



The stationary Ornstein-Uhlenbeck sheet {X(s,t) : s,t G 1R} defined by covariance 
structure (II. ip can be represented as 



X(s, t) = e - as "^ W (e 2as , e m ) , s,teR. 

Consider the sheet Z(s, t) :— mihi(s, t) + ■ — \-m p h p (s, t) +X(s, t), (s, t) G 
Theorem 12.11 for the functions 



, . 2y/a/3uv flogu logf 

9M = ^—h k 



k = 1,2,. ..,p, 



(2.4) 
Applying 

(2.5) 
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and for the domain G bounded by the curve T := r lj2 ur 1 ur 2 U r , where 

fi, 2 := {(«,7i») : u G [e 2aa ,e 2 ^]}, f x := {(u,7i(u)) : u G [e 2a \e 2QC ]}, 
T 2 :={(uMu)):ue [e^.e**]}, f := { {u, 70W) : u G [e 2 ^, e 2ac ] }, 

with 

7» :=exp (2/3 7j (log( U )/2a)), j G {{1, 2}, {1}, {2}, {0}}, 
we obtain the following result. 

Theorem 2.3 If hi, . . . , h p are twice continuously differentiable inside G and the partial 
derivatives d\hi, d 2 hi and d\d 2 hi, i = 1, . . . ,p, can be continuously extended to G then 
the probability measures Pz and P^, generated on C{G) by the sheets Z and X, 
respectively, are equivalent and the Radon- Nikodym derivative of Pz with respect to P^ 
equals 

dP z 



(Z) 



= exp I (m T Am - 2C T m) j , 



where A := (A M )^ =1 , m := (m 1 , . . . , m p ) T and ( := (0, . . . , C P ) T with 

A k ,e :=h k (a,'j 2 (a))h e (a, 72(a)) + h k (c, 71(c) )h i (c, 71(c)) 

+ fyfc(&i,7i(&i))M&i>7i(&i)) + M 02 ' 72(02)) h e (b 2 ,-f 2 (b 2 )) 
61 

+ y a/i fe (s,7i i2 (s))^( s ' 7i,2(s)) + a _1 9i/i fe (s,7i i2 (s))(9i^(s ) 7i,2(s)) 

a 

c 

+ y a/i fc (s,7i(s)) +9i/i fe (s,7i(s)) /^(s,7i(s)) + a _1 9i^(s, 7 i(s)) 



61 



&2 

+ y a/i fc (s,7 2 (s)) - 9i/ifc(s,7 2 (s)) /^(s, 72 (s)) - aT^/i^s, 7 2 (s)) 



ds 



ds 



ds 



(2.6) 



+ 



a/n(i,7o(i))ft<(s,7o(s)) +« 1 dih k (s, n/o(s))dih e (s, 7o(s)) 



ds 



62 



71,2(1) 



+ 



fih h (^(t),t) - d 2 h k (^l(t),t) h e fcl{t),t) - r 1 d 2 h e {%l(t),t) 



dt 



71,2(61) 
71(c) 



+ 



^ fe ( 7 r i (0,t)^(7r 1 (0^)+r 1 ^(7r 1 (0^)^(7r 1 W^) 



dt 



71 (61) 
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72 (62 ) 



+ 



/3h k (^\t),t)h e (^\t),t) + f3- l d2h k {^ 1 {t),t)d 2 h l {^ 1 {t),t) 



dt 



72(a) 
70(62) 



+ 



Ph k (% l {t),t)+d 2 h k (%\t),t) h t (%\t),t) +p- 1 d 2 h e ( 1 ^ 1 (t),t) 



dt 



70(c) 



+ // [a/3h k (s,t)h e (s,t) + a 1 (3dih k (s,t)dih l {s,t) + a(3 1 d 2 h k (s,t)d 2 h e (s,t) 



+ a (3 did 2 h k (s,t)did 2 h f (s,t)] ds dt 



G 



and 



C*:=4/ifc(&i,7i(&i))£(&i,7i(&i)) 



(2.7) 



+ 2 



a/i fc (s,7i(s)) + dih k (s^i{s)) Z(s, 7i(s))ds + a 1 Z(ds,j 1 (s)) 



+ 2 



a/ijfc(s,7i, 2 (s)) - <9ifyfc(s, 71,2(3)) Z(s,7i )2 (s))ds - a x Z(ds, 71,2(5)) 



72(62) 



+ 2 



72(a) 
71.2(a) 



+ 2 



Ph k (^(t),t) +d 2 h k ( 1 ^(t),t) Z( 7l - 2 1 (t),t)dt + /3- 1 ^(7i: 2 1 W,dt) 



+ / / [a/3/i fc (s,t) + /3dih k (s,t) + ad 2 h k (s,t) + d]_d 2 h k {s,t)\ 



71,2(61) 



G 



X 



[Z(s, £)ds dt + a _1 Z(ds, £)dt + 0~ x Z{s, dt)ds + a -1 /3 -1 ^(ds, dt)] . 



TTie maximum likelihood estimator of the parameter vector m based on the observations 
{Z(s,t) : (s,t) G G} has the form m = v4 -1 £ and has a p- dimensional normal distribution 
with mean m ano? covariance matrix A" 1 . 

Finally, consider the sheet Z(s, t) := mihi(s, t) + ■ ■ • + m p h p (s, t) + X(s, t), (s, t) G IR+, 
where {X(s,t) : s,t > 0} is the zero start Ornstein-Uhlenbeck sheet defined by ( 11. 2(1 . For 
a / and /3 7^ the sheet {X(s, t) : s,t > 0} can be characterized as a zero mean 
Gaussian process with 

•2 



EX{s 1 ,t 1 )X{s 2 ,t 2 ) = ^-(e 



). 
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hence, for example, in case a > and ft > it can also be represented as 



X(s,t) 



a 



2y/a~ft 



(2- 



In this way similarly to the stationary case one can apply Theorem 12.11 for the functions 

k = l,2,...,p, (2.9) 



, . 2y/aft{u+l){v + l) , f\og(u + l) log(v + l) 
#fc(w,t>) = — h k 1 



a \ 2a 2ft 

and for the domain G bounded by the curve Y := r 12 ur 1 ur 2 U r , where 
ri, 2 :={(n,7i, 2 H) -ue [e 2Qa -l,e 2a6l -l]}, fx := {(«, 7i(«)) : u G [e 2afel -1, e 2Qt - 1] } 



r 2 

with 



:={(w,7 2 (w)) :ue[( 



2 «a 



l,e 



2 «6 2 



1 



2ab 2 j e 2aC 1 



%(«):= exp(^2/3 7j (log( M + l)/2 
and obtain the following result. 



"o :={(w,7o(«)) : w e [ 
1, jg{{1,2},{1},{2},{0}} 



Theorem 2.4 If a ^ and ft ^ 0, functions hi,...,h p are twice continuously 
differentiable inside G and the partial derivatives d\hi, d 2 hi and d\d 2 hi, i — l,...,p, 
can be continuously extended to G then the probability measures Pz and Px, generated 
on C(G) by the sheets Z and X, respectively, are equivalent and the Radon- Nikodym 
derivative of Pz with respect to Px equals 



dP; 



-(Z) = exp (-^4(m T Am - 2C T m) 1 
[ 2a 1 J 



&P X 1 2(T 2 

where A := (A k/ ) P ke=1 , m := (m 1) . . . , m p ) T and ( := (Ci, ■■■,(?) ' with 

A k)l := coth (ad) coth (ftj 2 (a))h k (a, j 2 (a))h e (a, 72(a)) + h k (c, 71(c)) hi(c, 71(c)) (2.10) 
+ coth (/37i(6i))/i fc (6i,7i(6i))/i^(&i, 71(61)) + coth (ab 2 )h k (b 2 , 72 (6 2 ))/i£(6 2 , 72(62)) 

+ /coth (/3 7l)2 (s)) a/ifc(s, 7 i i2 (s))^(s, 7 i i2 (s))+a" 1 <9 1 /i fc (s, 7li2 (s))9 1 ^(s, 7li2 (s)) ds 



+ /coth (/3 7l (s)) a/ifc(s,7i(s))+9i/ifc(s,7i(s)) ^(s, 7 i(s))+a 1 d 1 h e (s, 7i(s)) ds 



61 



6 2 



+ 



a coth (as)/i fc (s,7 2 (s)) - dxh k (s, j 2 (s)) 



x 



coth (as)he(s, 72 (s)) — a 1 9i^(s, 72 (s)) 



d.s 
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ah k (s,~f (s))h e (s,7o(s)) +a 1 d 1 h k (s,-f (s))d 1 h e (s,-f (s)) 



ds 



62 



71.2(a) 

+ J coth (a^(t)) [/3coth (Pt)h k fcl(t),t) - d 2 h k (^(t),t) 

71,2(61) 



x 



coth {/3t)h e fal(t),t) - /3- 1 d 2 h e (^(t),t) 



dt 



71(c) 



+ 



Ph^^^^h^it),^ + f3- 1 d 2 h k (^ 1 (t),t)d 2 h e {^ 1 (t),t) 



dt 



71(61) 
72(62) 



+ / coth {a l2 \t)) (3h k ( l2 \t),t)h e ( l2 \t),t)+r 1 d 2 h k (^ 1 (t),t)d 2 h e ( l2 1 (t),t) 



dt 



72(a) 
70(62) 



Ph k ( lo 1 (t),t)+d 2 h k (% 1 (t),t)\ [h e (^\t),t) +p- 1 d 2 h e (^\t),t) 

+ J J [a(3h k (s, t)h e (s, t) + a^pfyh^s, t)dih e (s, t) + a/3~ 1 d 2 h k (s, t)d 2 h e (s, t) 
+ a^p^d^h^s, tfadzhtis, t)] ds dt 



dt 



70(c) 



G 



and 



C*:= 



+ 



61 



x 



61 



+ 



+ 



1 + coth(a&i)] [l + coth (/3 7 i(&i))] /ifc(6i,7i(6i))^(6i, 7i(6i)) 
1 + coth (/?7i( s )) «/?-fc(s,7i( s )) + d x h k {s, 71 (s)) 

Z(s,7i(s))ds + a _1 Z(ds,7i(s)) 
1 + coth (^7i, 2 (s)) a coth(as)/i fe (s, 7i, 2 (s)) - <9i/i fe (s, 71,2(5)) 
coth(as)Z(s,7i i2 (s))ds - a^Z(ds, 71,2(5)) 
1 + coth (a72 _1 (0)] [Ph k ( l2 -\t),t) + d 2 h k ( l2 \t),t) 
Z( l2 \t),t)dt + ^ 1 Z( l2 \t),dt) 



(2.11) 



x 



72(62) 



/ 

72(a) 



X 
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Figure 2: Means of the estimates of the components of m in Example 13 . 1 1 for 25 < n < 100. 



7i, 2(a) 



+ 



1 + coth (ar,$(t)) Ph k fcl(t),t) + d 2 h k { 1 ^ 2 {t),t) 



71,2(61) 



+ 



J j [a(3h k (s, t) + Pdxhkis, t) + ad 2 h k (s, t) + did 2 h k (s, t 



G 



[Z{s, t)ds dt + a~ 1 Z{ds, t)dt + 0~ x Z{s, dt)ds + a^/T^ds, dt)] . 



The maximum likelihood estimator of the parameter vector m based on the observations 
{Z(s,t) : (s,t) G G} has the form rh = A~ 1 ( and has a p- dimensional normal distribution 
with mean m and covariance matrix A^ 1 . 



Remar k 2.5 Observe tha t Theorems 12. 1[ [2731 and I2~4l are generalizations of Theorems 4, 5 
and 6 of lBaran et al\ ( 120031 ). respectively, where G — [Si, S 2 ] x [Ti , T 2 ] , with [Si , S?], [Ti , T o] C 
(0, 00). Hence, from these theorems one can also derive the results of lArato. N. M.I ( 1l997al lbl). 
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3 Simulation results 

To illustrate our theoretical results we performed computer simulations using Matlab 2010a. 
In order to simulate the Gaussian random fields considered above their Karhunen-Loeve 
expansions are applied. For the Wiener sheet W(s,t), < s < S, <t <T we have 

, , -A 8v^T /tt(2j-1)A /tt(2A;-1)s\ , , 



wher e {u}j,k '■ 1 < j,k <n} are independent standard normal random variables (IDeheuvels et al. 



200 



(]). The expansions for stationary and zero start Ornstein-Uhlenbeck sheets X(s,t) and 



X(s,t) with < s < S, <t <T can directly be derived from ( 13.1 p using representations 



f)2.4p and (12.81) . respectively (see e.g. IJaimez and Bonnetl (119871 )). yielding 



" 4ae a(s-s)+p{T-t) / ( 2 j _i) e 2^-r)\ /7r(2A;-l)e 2Q ( s - s )\ 

A 4a v /(e 2 " g -l)(e 2 ^-l) . f n(2j - l)(e 2/3t - 1) \ 

~ ~ Uhk ^^^{2k - l)(2j - 1) Sm V 2(e 2 ^ - 1) ) 

'7r(2fc- l)(e 2as - 1) 



x sm 



2(e 2a5 - 1) 
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Figure 4: Means of the estimates of the components of m in Example I3.2l for 25 < n < 100. 



In each of the following examples 1000 independent samples of the driving Gaussian 
sheet were simulated with n varying between 25 and 100 and the means of the estimates of 
the parameter vector m and the empirical covariance matrices of the vectors ( defined 
by (12. 3p . (12. 7p and (12. lip , respectively, were calculated. 



Example 3.1 Consider the model 



Z(s, t) = mi (s 2 + t 2 ) + m 2 (s + t)+ m 3 (s ■ t) + W(s, t), 



[s,t)eG, 



where W(s,t), (s,t) G [0,8] 2 , is a standard Wiener sheet and G is a circle with center 
at (6, 6) and radius r = 2. In this case the entries of the matrix A defined by (12. 2p 
and the approximations of the components of ( = (Ci ? C2, Cs) T defined by (12. 3p can be 
calculated using numerical integration, where Matlab function quad is applied (recursive 
adaptive Simpson quadrature). 



The theoretical parameter values are mi 
covariance matrix of ( equals 



5, m 2 = 8 and 777,3 = 3, while the theoretical 



.4 



/339.0895 38.6688 128.0000\ 
38.6688 5.9115 16.0000 
128.0000 16.0000 50.5752 
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Figure 5: Estimated covariances of ( in Example I3.2I for 25 < n < 100. 



On Figure [2] the means of the estimates of the three parameters, while on Figure |3] the 
estimated covariances of ( are plotted versus the rate n of the approximation ( 13. ip . In 
case of n = 100 we have (5.0157, 7.9868, 2.9594) for the mean and 



/338. 



A 



7473 38.8697 127.578 



4\ 



38.8697 5.9995 16.0409 
^127.5784 16.0409 50.2610 



for the covariance matrix. 



Example 3.2 Consider the model 



Z(s, t) = mi(s 2 + t 2 ) + m 2 (s + t) + m 3 (s ■ t) + X(s, t), 



{s,t)EG, 



where X(s,t), (s,t) e [0,3] 2 , is a stationary Ornstein-Uhlenbeck sheet with parameters 
a — 1, (3 = 1, cr = l and G is a circle with center at (2,2) and radius r = 1. Similarly 
to Example 13.11 the entries of the matrix A defined by ( 12 .6p and the approximations of 
the components of ( = (Ci)C2,C3) T defined by (12. 7D can be calculat ed using numerical 
integration, where Matlab functions quad and quad2d ( IShampind . 120081 ) are applied. 

The theoretical parameter values are the same as before, that is mi = 5, m 2 = 8 and 
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Figure 6: Means of the estimates of the components of m in Example I3.3l for 25 < n < 100. 

m 3 = 3, while the theoretical covariance matrix of ( equals 

/l797.8554 682.2301 801.6460^ 
A = 682.2301 274.1195 305.9734 . 
^ 801.6460 305.9734 382.9247y 

On Figure H] the means of the estimates of the three parameters, while on Figure the esti- 
mated covariances of ( are plotted versus the rate n of the Karhunen-Loeve approximation. 
In case of n = 100 we have (4.9938, 8.0098, 3.0030) for the mean and 

/ 1806.4240 681.5420 809.0529\ 

A= 681.5420 271.5779 304.4110 

^ 809.0529 304.4110 380.8634y 

for the covariance matrix. 



Example 3.3 Consider the same regression 

Z(s, t) = mi(s 2 + t 2 ) + m 2 (s + t) + m 3 (s • t) + X(s, t), 



(s,t)eG, 



as in Examples 13.11 and 13.21 but now the driving process X(s,t), (s,t) G [0,3] 2 , is a zero 
start Ornstein-Uhlenbeck sheet with parameters a — 1, /3 = 1, a — 1. Similarly to Example 
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Figure 7: Estimated covariances of ( in Example I3.3I for 25 < n < 100. 



13.21 G is a circle with center at (2, 2) and radius r — I. Again, the entries of the matrix 
A defined by (12.101) and the approximations of the components of ( = (Ci, (2, Cs) T defined 
by (12. lip can only be calculated with the help of numerical integration. 

The theoretical parameter values are the same as before, that is rri\ = 5, m 2 = 8 and 
m 3 = 3, while the theoretical covariance matrix of ( equals 



/l892.7035 725.4822 843.230l\ 



A 



725.4822 295.8952 321.8680 
843.2301 321.8680 395.0477y 



Similarly to the previous examples, on Figure |6] the means of the estimates of the three 
parameters, while on Figure [7] the estimated covariances of ( are plotted versus the rate n 
of the Karhunen-Loeve approximation. In case of n = 100 we have (5.0067, 7.9992, 2.9866) 
for the mean and 



.4 



/1884.8293 719.6687 845.4606\ 
719.6687 294.0491 321.0663 
840.0300 321.0663 396.9913 y 



for the covariance matrix. 
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A Appendix 

A. 1 Proof of Theorem 12^1 

To prove Theorem 12.31 one has to use representation (12. Ah and apply Theorem 12.11 for the 
random field 



Yin, v) := ^ = migi{% v) + _ + mMUt v) + w{% v) 



observed on G, where functions g k are defined by (12. 5p . As 



dig k {u,v) = — 7=(a + di)h k ' 



ay/au V 2a ' 26 



$,<&(«,«) = ^=(/3 + ft)/** (-2^, J, (A.l) 

^ / a \ (a /logM logf 



did 2 q k (u,v) = (a + d\)(B + d 2 )h k ( 

y V ' ; 2a v / ^m^ V 2a ' 2/3 



short calculations show 



A k ,e =4/ife(6i,7 li2 (&i)) /i£ (61,71,2(^1)) 



61 



(a- <9i)/i fe (s,7i )2 (s)) (1-a 1 9l)^(s,7i, 2 (s)) 



ds 



a + <9i)/i fe (s,7i(s)) (1 + a 1 d t )/i/(s,7i(s)) 



ds 



61 

72(62) 



+ 2 



(/3 + 9 2 )/i,(7 2 - 1 (t),0 (l + r 1 ft)^(7 2 " 1 (*) ) *) 



dt 



72(a) 
71,2(1) 



+ 2 



(/3 + d 2 )h k (j^(t),t) (1 + /3- 1 9 2 )/ i ,( 7l - 2 1 (t),t) 



dt 



71,2(61) 



+ 



a + d x ){B + d 2 )h k (s, t) (1 + a- 1 d l )(l + d-^htis, t) 



ds dt 



G 



-4} + aA% + /Mg + a-USJ + /T^g + Ag, 
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where 



A$ :=Ah k (b 1 niMhe(biniM)+ 2 J d ± [^(s,7i(s))^(s,7i(s)) 

b 1 72(62) 



ds 



2 J d 1 [h k (s ni , 2 (s))h e (s ni , 2 (s))\ds + 2 J d 2 [h k (^ 1 (t),t)h e (^ l (t),t) 

72 (a 



dt 



71.2(a) 



+ 2 



J d 2 [h k ( 1 ^(t),t)hz( 1 ^(t),t)~\ dt + J J ^[hkis^h^t)] dsdt, 

9.(h) G 



71,2(61) 



A ki'-= 2 J h k (s,-f 1 (s))h e (s,-f 1 (s)) ds + 2 J h k (s,j lt2 (s))h e (s,-f h2 (s)) ds 

bi a 

+ JJ d 2 [h k (s,t)h e (s,t)] dsdt, 



G 

72(62) 



71.2(a) 



4 3 j:=2 J h k (^ 1 (t),t)h £ (~f 2 1 (t),t)dt + 2 J h k (^(t),t)h £ (^l(t),t)dt 

72(a) 71,2(61) 

+ JJ 3i[/i fc (s,f)MM)] dsdt, 

G 

c 61 

Aj$:=2 y 9i/i fe (s,7i(s)) 9i^(s,7i(s)) ds + 2 ^ <9i/i fe (s, 71,2(3)) dxh^s, 7i, 2 (s)) ds 

61 a 

+ JJ d 2 [d 1 h k (s,t)d 1 h e (s,t)] dsdt, 



G 
72(61) 



71.2(a) 



2 y 9 2 ^.(72 _1 W^)«92^(72 _1 W ) dt + 2 y d 2 h k (^{t),t)d 2 h i (^l(t),t)dt 

72(a) 71,2(61) 

+ ^ d 1 [d 2 h k (s,t)d 2 h e (s,t)] dsdt, 

G 

Afl := JJ [a(3h k (s, t)h e (s, t) + a^^h^s, t^h^s, t) + a^d^s, t)d 2 h e (s, t) 



+ a^p^d^hkis, ijdidzheis, t)] ds dt. 
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Obviously 



b 2 &i 

d 2 [h k (s,t) h e (s,t)]ds dt = / /i fc (s,7 2 (s))/i < (s,7 2 (s)) ds - / /i fc (s, 7i, 2 (s))/^(s, 71,2(5)) ds 



-y ^fc(s>7i(*))^(s,Ti( s )) ds + J h k ( s ,lo(s))h e (s,j (s)) ds, (A. 2) 

61 b 2 

71(c) 71,2(1) 

^LpifcCM) h t {s,t)]dsdt = J h h (ji\t),t)hi(j{\t),t) dt -J h k (^t\t)h t {^{t\t) dt 

G 11(h) 71,2(61) 

72(62) 70(62) 

- J h k (72 -1 (t) , t) ft/ (7 2 _1 (*) ,t)dt + J h k (7^ (t) , t) ^ ^ (t) , t) dt , ( A.3) 

72(a) 70(c) 

and similar expressions can be derived for 

J J d 2 [dih k (s,t) dihi(s,t)] ds dt and J J d%[d 2 hk{s,t) <9 2 /^(s,t)] dsdt, 



respectively. Hence, using also that 
2 ^ [/i fe (s, t)h e (s, t)] ds dt 



(A.4) 



G 



l>2 



3| 



72(62) 

ds- / a 



^* (7a _1 (*). *)^/ (72 1 



dt 



72(a) 



71 (c) 



-y* 

61 
61 

- /ft 



/i fc (s,7i( s ))^( s ; 7i(s)) 



^*(sj7i,a(s))^(s,7i )a (s)) 



ds+ y <9 2 

71(61) 

71,2(61) 

ds + J 3 2 

71.2(a) 
70(c) 



/ ifc ( 7 r 1 (t),t)^( 7 r 1 (t),t) 



dt 



dt 



+ y Sl 

62 



/ifc(s,7o(s))^( s '7o(s)) 



ds- y <9 2 

70(62) 



/i fc ( 7o - 1 (t),t)^( 7o - 1 (t),t) 



dt, 



after tedious but straightforward calculations we obtain (12. 6ft . 
Further, using again (lA.ljl and Remark 12.21 we have 



t - /•(!) 1 /-(a) 1 /-(a) 1 X 4 ) 1 /^( 5 ) 1 X 6 ) 

SA: — Sfc ~T Sfc "I" Sfc < S>k "1" V>fc "r Sfc 5 
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where 



A (2) 



= 4/i fc (6i,7i^(6i))Z(6i, 71,2(61)) 

c 

l.i.m.- / e as (a + d 1 )h k (s, ll (s)) 
e->o gj 

■ (e „, + g) -A z ( M^ri 7i(8) ) _ <- z( , i7lW) 



ds, 



a(3) 



4 y (a - 9i)/i fc (s,7i )2 (s))Z(s,7i i2 (s))ds 

a 

61 

-l.i.m.- / e as (a- ^)^(s, 7l)2 ( S )) 
e^o Q J 

a 

■ (e2 «, + e) ^( l0g(e2 2 " +g) , 7l , 2 ( S )) - e°»Z( S , 71 , 2 ( S) ) 



ds, 



72(^2) 



Sfe 



l.i.m. - 

5^0 5 



(P + d 2 )h k ( l2 1 (t),t) 



72(a) 



lu . log(e 2 ^ + 5) 



2/3 



dt, 



71.2(a) 



A (5) 



l.i.m. - 

5^0 (5 



B /3 '(/3 + 5 2 )/ ifc (7i: 2 1 W,t) 



71,2(^1) 



X 



(e^ + 6) 1/2 z[ % l(t), 



lu . log(e 2 ^ + 5) 



AG) 



l.i.m. — lie 

g,s^o go 

G 



as+f > t ^ + d 1 )(P + d 2 )h k (s,t) 



la 



2/3 

e as+ ^Z(s,t) 



^z(^l(t),t) 



dt, 



{e 2as + Q ) l,2 ^ t z( K 



2a 



t)-e as e 



+ (e 2as + Q) 1/2 (e 2 ? t + 5) 1/2 z( 



\og(e 2as + q) log(e 2/3 * + 6) 



2a 



2/3 



ds dt, 



implying ( 12. 7p . 
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A. 2 Proof of Theorem I2T41 



The proof, which we give only in the case a > 0, > 0, is similar to the proof of Theorem 
12.31 Here one has to use representation (12. 81) and apply Theorem 12.11 for the random field 



Y(u,v) 



2y/aP(u + l)(v + l) Aog(u + l) \og(v + 1) 

a V 2a ' 2/3 

■ rrngtiu, v) + . . . + m p g p (u, v) + W(u, v) 



observed on G, and functions g k should be defined by ( I2.9p . In this way 

dig k (u,v) 



+ 1) / , o U fl°g(u + l) log(« + l) 



d 2 g k {u,v) 
did 2 g k {u,v) 



oyJa(u + 1) 

+ 1) 



1 



2a ' 2/3 
log(w + 1) log(f + 1) 



2a 



2/3 



(A.5) 



2ay/a/3(u + l)(v + l) 



(a + d 1 )(f3 + d 2 )h k 



\og(u + 1) log(v + 1) 



2a 



2/3 



so 



.4 



k.C 



e 2a&! e 2/3 7 i, 2 (fei) 



l e^7i, 2 (6i)_i fefc ( 6 i' 71,2(61)) M fe i> 71,2(61)) 



bi 



+ 2 



+ 2 



+ 2 



3 2/37i,2 (a) 
e 2/3 7 i, 2 (s)_l 

e 2/37i0) 

e 2/3 7 i(s)_l 
&l 

72(62) , 

e 2o 72 *(t) 



(acoth(as) — 9i)/i fe (s,7i )2 (s)) (coth(as)-a 1 9i)/i < (s, 71,2(5)) 



ds 



a + 9i)/i fe (s,7i(s)) (1 + a 1 9 1 )/i < (s,7 1 (s)) 



ds 



72(a) 
71,2(0) 



e 2a 72 "" 1 (t)-l 



(/3 + d 2 )h k {^ 1 (t),t)} Ui + p-^M^it),^ 



dt 



+ 2 



71,2(61) 



e 2«7i"2(*) 
3 2a 7 - 2 1 (t)_ 1 



(/3 + 5 2 )^( Tl ; 2 1 (t),t)l [(1 + 0-%M%l(t),t) 



at 



+ 



(a + 0009 + d 2 )h k {s, t)} \{1 + a~ 1 0O(l + /3" 1 2 )^(s, t) 



ds dt. 



G 



Collecting separately the terms containing 7 12 , 71, 7 2 and 70, after long straightforward 
calculations using again (I A . 2 1) f |A~4|) we obtain 



Al t _ + 4(2) , A (3) , 4(4) , A (5) 
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where 

A k\ '■= \ h k {K 71,2(61)) hi (b u 71,2(61)) ~\ h k (a, 71.2(a)) ^ (a, 71.2(a)) 
+ coth (era) coth (/3'fi :2 (a)) h k (a, "fi j2 (a)) hi(a, "fi :2 (a)) 

+ coth (^71,2(60)^(61, 7i,2(6i))/i<(6i, 71,2(61)) + coth (aa)h k (a, ~f lj2 (a))h e (a, 71,2(0)) 
61 

+ / coth (/37i, 2 (s)) aa fc (s,7i i2 (s))a <! (s,7i i 2(s)) + a _1 9i/ifc(s,7i i2 (s))9i^(s ) 7i,2(s)) ds 



71.2(a) 

+ J coth (cry^*)) [(/?coth(/3t) - d 2 )h k (^{t),t) 



l(4) 



71,2(^1) 



4 (2) ._ 



(coth(^) - /T 1 ^) a, fe 1 (*),*) 
j , 1 

-a fc (c,7i(c))/i £ (c, 71(c)) - -/ife(6i,7i(6i))^(6i,7i(6i)) 



df, 



+ 



y coth (/?7i(s)) (a + c\)/i fc (s,7i(s)) (l + aT 1 ^)/^, 7i(«)) ds 



71 (c) 



+ 



Ph k (^\t),t)h t (^\t),t) + f3~ 1 d 2 h k (^ 1 (t),t)d 2 h e (^\t),t) dt, 



71(61) 



^Lj ■= ^h k {b 2 ^ 2 {b 2 ))h t {b 2 , 72(62)) - -h k (a,j 2 (a))h e (a, 72(a)) 

+ coth (062)^(62,72(62)) ^(62,72(62)) - coth (aa)/i fc (a,7 2 (a))/^(a,7 2 (a)) 



62 

+ J (acoth(as) - d 1 )h k (s,j 2 (s)) ( coth(as) - a' 1 d 1 )h e (s, 72(5)) 



ds 



a 

72(62) 



+y"coth (a 72 - 1 (t))[/3/ ifc (7 2 - 1 (t),t)^(7 2 - 1 (t),t) + p- 1 d 2 h k ( 1 ; 1 (t),t)d 2 h £ ( l2 - 1 (t),t)]dt, 

72(a) 

1 1 

= o^( c '7o(c))a^(c, 70(c)) - -h k (b 2 , 70(62)) ^(62,70(62)) 



ds 



+ y a/ifc(s,7o(s))/i/(s,7o(s)) + a 1 d 1 h k (s, 7o(s))<9A( s > 7o(s)) 
62 

70(62) 

+ y [(/3 + 92)a fe (7o~ 1 W^)][(l + r 1 92)a,(7 - 1 (t),t)]dt ) 

70 (c) 

4 5 J := y^ [a/9/i fc (s, t)a,(s, t) + a^^h^s, tfahfa, t) + a/T^/i^s, t)d 2 h e (s, t) 

G 

+ a~ 1 /3~ 1 (9i(92a fc (s, t)did 2 he(s, t)] ds dt, 
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which equals (I2.10p . 

At the end, using again (1A.5j) and Remark 12.21 we have 



t - A 1 ) A 2 ) , A® _j_ a(4) , AS) >(6) 
SA: — Sfc i Sfc "I" Sfc "r Sfc "r Sfc *r" Sjfc > 



where 



Sfc 



a(2) 



A (3) 
S>fc 



a(4) 
Sfc 



4 e 2a6i+2 i a 71 , 2 (bi) 



e 2abi _ lV e 2^ 7l , 2 (6i) _ 1) 
c 

1 r 4e as+2/371 ( s ) 



/ife(6i, 7i,2(6i))^(&i, Ti,a( 6 i)) 



l.i.m 

e^o g J e 

61 



.2^W-l (Qi + gl) ^( g ' Tl(g) ) 

(e^ + g )^( j^ , 7l (,))-e^( S , Tl (,)) 



ds. 



6l 



4 e 2as+2/3 7 i, 2 (s) 



(e 2QS - l)(e 2 ^i l2 ( s ) _ i) 



(a coth(as) - di)h k (s,j 1>2 {s))Z(s,j 1>2 {s))ds 



l.i.m. - 

£>->o g J e 



1 /" 4e as+2 ^' yi ' 2 ( s ) 



,2/3 7 i, 2 (s) 



- (a coth(as) - di)h k (s, 7i,2(s)) 



x 



2a xi/2„/log(e 2as + q) \ a . s 

6 Q > V 2^ >7i,2(s)J-e Z(s,7i i2 (s)j 



ds, 



72(^2) . 
1 /• 4 P /?*+2a72~ (*) 

'■Ho- ? J 5=?«-r&' + '(*)■«) 

72(a) 



(e^ + 5) 1/2 Z( 72 - 1 (t) 



log(e 2 ^ + <5) 



(5) ._ 



2/3 



71,2(1) 

1 /• 4e' 3 *+ 2a Ti,2(*) 
l.i.m. - / — (/9 + 9 2 )^(7i;2 1 (*),*) 



^(7^(0, t) 



dt, 



71,2(61) 



log(e 2 ^ + (5) 



2/3 



3^(7^),*) 



dt, 
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Cf := Li.rn.-i // e^ t (a + d 1 )^ + d 2 )h k (s,t) 



<5->0 q8 



G 



- (e- + g) 1/2 e^ (M^_±A, t ) _ e «- {e W +6) V* Z ( s , lggg±g 



+ (e 2as + f? ) 1/2 (e 2 ^ + 5) 1/2 z( 



2/3 

2^/log(e 2as + g) log(e 2/3t + <*)• 



ds dt, 



2a 2/3 

implying (EOT]) . □ 
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